Large deviations and heterogeneities in a driven kinetically con- 
strained model 
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Abstract. - Kinetically Constrained Models (KCMs) have been widely studied in the context of 
glassy dynamics, focusing on the influence of dynamical constraints on the slowing down of the 
dynamics of a macroscopic system. In these models, it has been shown using the thermodynamic 
formalism for histories, that there is a coexistence between an active and an inactive phase. This 
coexistence can be described by a first-order transition, and a related discontinuity in the derivative 
of the large deviation function for the activity. We show that adding a driving field to a KCM 
model does not destroy this first-order transition for the activity. Moreover, a singularity is also 
found in the large deviation function of the current at large fields. We relate for the first time this 
property to microscopic structures, in particular the heterogeneous, intermittent dynamics of the 
particles, transient shear-banding and blocking walls. We describe both the shear-thinning and 
the shear-thickening regimes, and find that the behaviour of the current is well reproduced by a 
simple model. 



Introduction. — Disordered materials, such as gels, 
glasses, complex fluids or granular materials, exhibit slow 
relaxation and heterogeneity, both in space and time. In 
particular, spatial dynamical heterogeneities, i.e the coex- 
. . , istence of mobile particles and blocked particles, have been 
^ ' widely studied jJL, and have benefited from the results of 
k>( KCMs (Kinetically Constrained Models), which by defi- 
5^ nition rely on specific dynamical rules depending on the 
, mobility of the particles [5] . 

When driven by an external force, these complex sys- 
tems exhibit specific spatial features, such as shear- 
banding [3.:- Another complication arises [4] with shear- 
thickening, where viscosity increases when the driving 
force increases [S]. In this paper, we study a KCM driven 
by an external field, using both global observables (large 
deviation functions) and a microscopic description. This 
model, first introduced by Sellitto [6], has the advantage 
of displaying both shear-thinning and shear-thickening 
regimes: at small field, the current of particles is propor- 
tional to the field, whereas at large fields, the system has 
a negative differential resistance, i.e the current decreases. 
Similar effects have also been observed in electronic sys- 



tems [7] where negative differential conductivity can occur 
under strong electric fields. 

In this context we have studied the thermodynamics of 
histories in the approach developed by Ruelle [HUH] for 
deterministic dynamical systems. This formalism can be 
extended to Markov chains [TU] , and has been already de- 
veloped for KCMs [11]. In these latter studies, where the 
observable is the activity, one is able to probe trajectories 
in phase space with either finite, or almost zero activity, 
according to a "chaoticity temperature", s, which is in- 
troduced in the formalism. Contrarily to systems without 
kinetic constraints, there is a first-order transition at s = 
in the activity, which translates into a discontinuity of the 
derivative of the associated large deviation function. This 
transition is the global signature of the observation of dy- 
namical heterogeneity. In the case of driven KCMs, we 
will focus on the large-deviation functions for both the ac- 
tivity and the current. We will also present a microscopic 
analysis of the structures involved in the dynamics, and 
relate them to global quantities. 

Ruelle's formalism. — This formalism yields infor- 
mation about the fiuctuations of temporal trajectories 
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in configuration space, while the usual canonical ther- 
modynamics approach yields information about the fluc- 
tuations associated with configurations. One defines a 
dynamical partition function, with a chaoticity temper- 
ature that is conjugated to the observable. Hence, if O 
is a time-extensive observable, Zo{s,t) = ^e""*^). The 
brackets mean that one is averaging over all possible 
histories between and t. More precisely, Zo{s,t) = 
S/iistories P'"ob(history)e~^°('"*'*°''y\ where a history is 
the time series of configurations {Ci, C2, ■ ■ ■ , Ct}. In the 
long time limit, the partition function behaves like 

lim Zo{s,t) oc exp {i-0o(s)} 

where tpo is the large deviation function associated with 
the observable O. In practice, one is also interested by the 
mean of the observable in s-space, i.e, where the trajecto- 
ries are biased by the factor e""*^. This quantity is related 
to the large deviation function through the relation, in the 
large time limit: < 0{s,t) > /t = —iPq{s). 

The model. — In this paper, we study a two dimen- 
sional asymmetric simple exclusion process (ASEP) on a 
square lattice with dynamic constraints [6]. One can also 
view the model as the KA (Kob- Andersen) model in two 
dimensions -an example of KCM- |12j , driven by an exter- 
nal field. The dynamical rule is such that a particle can 
move to a neighboring site only if it has at least two empty 
neighbours, before and after the move. Such a constraint 
mimics the fact that a particle in a dense disordered ma- 
terial can only move if it has enough free volume. Such 
a model always has a non-vanishing diffusion coefficient 
except for a density equal to 1 [13]. The motion of the 
particles is biased in one direction of the lattice by an ex- 
ternal field E, such that the probability for a particle to 
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move in a certain direction is p = min 1 1 , exp 

where r indicates the displacement vector. The boundary 
conditions are periodic in both directions. In [5], the cur- 
rent was measured as a function of density p and the field, 
showing several distinct behaviours. For densities below 
a threshold value (approximately pc « 0.79) the system 
shows a monotonic growth of the current as a function 
of the external field (like in the non constrained ASEP 
model) , which eventually saturates to a maximum current 
value. For p > Pc, the current has two different regimes: 
for E small, a first ohmic regime (also referred to as a 
" shear thinning" regime) where the current grows propor- 
tionally to the field; and for E large, a negative differential 
resistance regime (also called "shear thickening" regime) 
where the current decreases as the field increases. Such 
a behaviour has also been observed in other similar mod- 
els of glasses [M]. In this paper, we are mainly interested 
in the origin of the non-monotonic behaviour of the cur- 
rent as a function of the field. As a consequence we will 
consider only densities above the critical density pc- 

In the following, we will consider two different time ex- 
tensive observables, namely the activity K and the in- 
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Fig. 1: The mean activity K{s) in the modified dynamics, nor- 
malised by its value in s = (top), and its large deviation 
function tpK{s) (bottom) for a system of size = 60^ at high 
density p = 0.80. For the simulation, we have considered 300 
clones of the system evolving for a total time of r = 5000 MCs. 
For any value of the external field there is a discontinuity in the 
derivative of ipK which corresponds to two different regimes in 
K{s)/t: an active one (s < 0) and an inactive one (s > 0). 



tegrated current Q. The activity K{t) is defined as the 
total number of moves from time to time t. The inte- 
grated current Q{t) — /p J{t')dt' is defined as the number 
of moves in the direction of the field from time to time 
t. The current calculated in 6^ is simply Q{t)/t. 

In order to simulate the biased dynamics (the "s- 
dynamics") and to compute the values of K{s), Q{s) and 
their large deviations functions tpxis) and i1^q{s), we have 
used the discrete time cloning algorithm proposed in ref- 
erence |15| . which was proved to be powerful enough for 
numerical studies of KCMs. 

Large deviations computations. — The large devi- 
ation function for the activity has already been studied for 
different KCMs without driving such as the Fredrickson- 
Andersen Model, the triangular lattice gas model, the 
Kob- Andersen model [TT]. A first-order phase transition 
in K was observed at s = 0. Such a transition is the 
global signature for the coexistence at s = (which corre- 
sponds to the unbiased, physical system) of histories that 
have positive values of the activity (typical of s < 0), 
and histories that have sub-extensive activity (typical of 
s > 0). It is a translation in phase space of the dynamical 
heterogeneities one finds when looking at configurations. 
For non-driven models the activity's phase transition is 
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Fig. 2: The mean integrated current Q{s) in the modified dy- 
namics normalized by the value Q{s = O)(top) and its large de- 
viation fimction iPq{s) (bottom) plotted with a proper rescal- 
ing. The system has a size of = 60^ sites and a density 
p = 0.8. The maximum current phase is attained for a field 
E* ~ 2. Here we have considered 300 clones of the system 
evolving for r — 5000 MCs. The increase of the discontinuity 
of Q{s) at s = is emphasized in the inset. 



directly linked to the presence of kinetic constraints, since 
it disappears as soon as the constraints are removed. In 
our driven model, we find that this first-order transition 
still exists for the activity (see Fig. [TJ, and is by no way 
smeared out by the driving, whether the system is in the 
positive or in the negative differential resistance regime. 
The discontinuity observed at s = increases as E gets 
larger. We checked that for all fields, the transition exists 
at the thermodynamic limit: we performed a finite-size 
study showing that the discontinuity is more and more 
visible as the size of the system increases. In order to do 
this study accurately, we needed more and more clones 
and longer simulation times for larger system sizes, in or- 
der for the cloning algorithm to converge. 

We considered then the integrated current, flowing in 
the direction of the external field. It is linked to the ac- 
tivity since it is the difference between the activity in the 
positive direction of the field and the one in the negative 
direction. In the case of a non-constrained asymmetric 
exclusion process in one dimension, it has been shown 
[16] that the current's large deviation function is a reg- 
ular convex function with a symmetry around E/2 due to 
the fluctuation theorem [T7]. This means that no discon- 
tinuity exists for any value of s for the derivative of ipg (s) 
in the non constrained case, or, in other terms, Q{s) is a 
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Fig. 3: The large deviation function ^q(s) for a system size of 
— 30^ sites and a density p = 0.8. Here we have considered 
300 clones evolving for different times r. We also show the re- 
sult for 1000 clones and a very long simulation time, suggesting 
the asymptotic symmetric shape, and a clear discontinuity of 
the derivative at s = 0. 



continuous function. We checked this in the case of the 
two dimensional ASEP. Moreover, a perfect collapse of the 
Q{s)/t and i^Q{s) curves can be found when rescaling s 
by s/ E and the ordinates according to the scaling used in 
Fig.m 

If we introduce the dynamical constraints we find (see 
Fig. corrections to this perfect scaling. For small fields 
no discontinuity can be found (even if one uses more clones 
and longer times in the simulation); the large deviation 
function V'q(s) has the same shape we should expect for 
a non constrained model (2d-ASEP). Conversely, if we 
increase the field, a gap between the left and the right 
derivatives of iPq{s) is generated, and the size of the gap 
increases continuously as we increase the field. In the limit 
of large system sizes, the discontinuity in the derivative of 
iPq{s) corresponds to a step discontinuity in the biased 
current Q{s); we also performed this finite-size study, and 
found that the discontinuity is sharper and sharper as the 
size increases. Here, the convergence of the cloning algo- 
rithm at fixed L is rather slow. Figure [3] illustrates this 
point: as time and the number of clones increase, the tran- 
sition is more and more visible. This figure also suggests 
that TpQis) should asymptotically converge to a symmet- 
ric shape with a central part equal to zero; this is actually 
what is predicted by the fluctuation theorem [17] and was 
found in a mean-field model |20j . This asymptotic curve 
could not be obtained here, due to numerical limitations: 
this difficulty was already raised by Sellitto [3] who also 
found deviations from the fluctuation theorem, even at 
very large times. 
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Fig. 4: On the left and middle: velocity field plots for two values of the external field, corresponding to the same value of the 
average current. The diffusive homogeneous behaviour at small fields can be contrasted with the localisation of transport in 
bands for large fields. On the right: the mean longitudinal velocity versus the direction transversal to the external field, for the 
same two values of the driving force. 



Microscopic analysis. — Physically, this transition 
in Q{s) and tpgis) at s = corresponds to a coexistence 
between histories that have a positive current (Q(s) is pos- 
itive for s < 0) from histories that have a subextensive 
current (for s — )> 0^ we have no current, and nearly all 
particles are blocked). 

In order to better investigate the emergence of the neg- 
ative differential resistance behaviour, we have analyzed 
some microscopic aspects of the dynamics. These micro- 
scopic observations are also useful in order to understand 
the behaviour of the large-deviation functions described 
above, in the regimes of small and large forcing. 

Simple observations of the dynamics of the system show 
that the driving by the field induces the formation of pe- 
culiar void structures, composed of quasi-linear segments 
of empty spaces, in the direction perpendicular to the 
field. The particles actually accumulate in some dense, 
blocked regions of the sample, separated by these "domain 
walls" . These domain walls are hard to remove because 
of the kinetic constraint; their lifetime is long compared 
to the microscopic dynamics. In other words, the driving, 
combined with the dynamical restrictions, is at the origin 
of structures that slow down the dynamics. This effect 
is enhanced when the field is increased, yielding shear- 
thickening behaviour. The domain walls can be directly 
visualized in the insets of figure[S]for small and large fields. 

At the level of single particle displacements, the exis- 
tence of such domains is also at the origin of different 
kinds of moves. For small fields (and small currents) the 
behaviour is almost diffusive, the particles randomly ex- 
plore the region around their starting point, even coming 
back and forth on their own path. In the negative resis- 
tance regime (large fields) few particles manage to cover 
very long distances with almost ballistic trajectories, while 
others remain caged in the vicinity of their starting point. 
The moving particles actually diffuse in the vertical di- 
rection for some time before making some direct jump in 



the direction of the field. Such jumps actually correspond 
to the sudden disappearance of the empty domain walls. 
This behaviour is probably at the origin of superdiffusiv- 
ity [B], which has also been observed in granular jammed 
systems [181. Moreover, it corresponds to a clear inter- 
mittent dynamics where the current is locally large for a 
short period of time. 

We have computed the distribution of wall sizes P{w) 
in the stationary regime. We have found that P{w) is well 
approximated by a simple exponential in a large range 
of values of p and E, allowing to define an average size 
< w > {p,E). This average wall size grows as the density 
or the field grows, and saturates at large values of the ex- 
ternal field [19] .We can in fact rationalize the behaviour of 
the current, J{E), by the following phenomenological ar- 
gument. Let us assume that one can attribute the slowing 
down of the dynamics by the mere presence of the walls. 
The average current counts the number of particles suscep- 
tible to move in the direction of the field, and not blocked 
by the presence of a wall: J{E) = A{l~e^^){l—pbiocked)- 
We assume that the probability to be blocked is propor- 
tional to the number of blocked particles confined on the 
left border of a wall. To lowest order in the wall size we 
assume a linear dependence of puocked on < w > (p,E): 
Pbiocked = a < w > {p, E), whcrc a depends in principle 
on p and E. We found that the curves for J{E) could be 
very well fitted by this formula, taking a independent of 
E (see Fig. [S]). The parameter a, which can be shown in 
a more refined study jT9] to depend on microscopic geo- 
metric details, such as the number of walls or the average 
distance between them, is found to be very weakly depen- 
dent on p in the range studied. The assumptions made 
here can explain the fact that the fit actually begins to 
fail at large p: then < w > gets larger, P{w) starts to 
deviate from an exponential form, showing that the ap- 
propriate expression for puocked may change. 

At the level of the entire sample, we have drawn the 
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velocity field plot for the small and large values of the 
driving force, from which one can obtain the velocity pro- 
file for the longitudinal component of the speed vector (see 
fig. 2]). Velocities in these plots are defined over a time 
lag of the order of the relaxation time; they also corre- 
spond to displacement snapshots on this timescale. These 
plots show that at large fields the system gets more and 
more organized, with bands of moving particles coexisting 
with static regions blocked by walls, reminiscent of shear 
banding observed in complex fiows. This illustrates the 
first-order transition for the large deviation of the current 
at large fields. The high velocity bands have a macro- 
scopic size and involve the whole longitudinal length of 
the system. At small field, the current is more homoge- 
neous, which is consistent with the absence of transition 
in the large-deviation function for the current. 




Fig. 5: Simulated flow curves (white squares) fitted by the 
phenomenological model J{E) = A{1 — e^^)(l — a < w >) 
(black dots). Densities range from 0.795 to 0.825, L = 100, 
and < w > is determined independently from the simulation. 
Below are shown two fractions of snapshots of configurations 
taken at p = 0.82, a) E=0, b) E=5, illustrating the presence of 
walls. The empty sites are in black, the particles in white: the 
empty sites organize into vertical walls at large fields. 



Conclusion. — We have found that the coexistence of 
mobile and blocked trajectories in configuration space is 
a conserved feature in KCMs even in the presence of ex- 
ternal driving. This is illustrated quantitatively through 
the study of the large deviations for the activity and for 
the current. In the shear-thinning regime, the microscopic 
behaviour is reminiscent of the dynamical heterogeneities 
found in KCMs, coming from dynamical restrictions; there 
is a transition in the large deviation function for the activ- 
ity, but not for the current. In the shear-thickening regime, 
the field induces a structuration of the flow, with domain 
walls separating dense jammed regions, as well as tran- 
sient shear bands. The discontinuity found at large fields 
in iPq{s = 0) illustrates the dominant role of the blocking 



walls in the structuration of the displacement field. We fi- 
nally note that we have very recently found a preprint by 
Speck and Garrahan on a related subject, where a first- 
order transition is also found for the entropy production 
for driven KCMs 
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